A first-principles model of time-dependent variations in transmission through a 

fluctuating scattering environment 
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Fading is the time-dependent variation in transmitted signal strength through a complex medium, 
due to interference or temporally evolving nmltipath scattering. In this paper we use random matrix 
theory (RMT) to establish a first-principles model for fading, including both universal and non- 
universal effects. This model provides a more general understanding of the most common statistical 
models (Rayleigh fading and Rice fading) and provides a detailed physical basis for their parameters. 
We also report experimental tests on two ray-chaotic microwave cavities. The results show that our 
RMT model agrees with the Rayleigh/Rice models in the high loss regime, but there are strong 
deviations in low-loss systems where the RMT approach describes the data well. 



Considering wave propagation between a source and 
a receiver in a complex medium, fading is the time- 
dependent variation in the received signal amplitude as 
the scattering environment changes and evolves [ij. Fad- 
ing is a challenging problem in many situations where 
waves propagate through a complicated scattering en- 
vironment. A common example is the nighttime varia- 
tion of AM radio signal reception in the presence of ray 
bounce(s) off a time varying ionosphere. Another com- 
mon observation of fading is experienced by radio lis- 
teners in automobiles moving among vehicles and build- 
ings in an urban environment. Fading exists in closed or 
open scattering systems and in all types of wave propaga- 
tion, and it is broadly studied in wireless communication, 
satellite-to-ground links, and time-dependent transport 
in mesoscopic conductors [ll-Ql- 

The fading amplitude is defined as the ratio of the re- 
ceived signal to the transmitted signal. The traditional 
models [l| of fading work well in certain regimes of ra- 
diowave propagation applications, where different prob- 
ability distribution functions are chosen depending upon 
the circumstances. However, these models are empiri- 
cally designed for particular scattering environments and 
frequency bands, and different (apparently unrelated) 
fitting parameters are introduced in different models. 
For example, the Rayleigh fading model applies a one- 
parameter Rayleigh distribution for the fading amplitude 
in an environment where there is no line-of-sight (LOS) 
path between the transmitter and the receiver, such as 
mobile wireless systems in a metropolitan area [l|-13| • The 
Rice fading model, on the other hand, applies a two- 
parameter distribution for situations with a strong LOS 
path [il, [a, Q . The detailed physical origins of these mod- 
els, and their parameters, are not clear. 

The complexity of the wave propagation environment 
is advantageous from the perspective of wave chaos the- 
ory because it means that wave propagation is very sensi- 
tive to details, and a statistical description is most appro- 
priate. For applying wave chaos approaches, the system 
should be in the semiclassical limit where the wavelength 
is much shorter than the typical size of the scattering sys- 
tem Q . Researchers have applied random matrix theory 



(RMT) in wireless communication [9| and analyzed the 
information capacity of fading channels |10l - [l3 |. or the 
scattering matrix (S) and the impedance matrix (Z) of 
the scattering system 13-15] . Here we directly apply the 
random matrix approach to the fading amplitude. 

We derive a RMT-based fading model that includes 
the Rayleigh and Rice fading models in the high loss 
regime, but the RMT model also works well in the limit 
of low propagation loss. In addition, the RMT approach 
combined with a model of non-universal features reveals 
the precise physical meanings of the fitting parameters 
in the Rayleigh and Rice models. 

Considering a 2 x 2 scattering matrix S which describes 
a linear relationship between the input and the output 
voltage waves on a network, the two ports can be assumed 
to correspond to the transmitter and the receiver. The 
fading amplitude is equivalent to the magnitude of the 
scattering matrix element 1521 1- We start with a RMT 
description of the 2x2 universal scattering matrix Srmt 
in a wave chaotic system [IJ]. This description assumes 
total ergodicity and does not account for system specific 
information (such as the coupling of the ports and the 
short ray paths between the ports). For time-reversal 
invariant (TRI) wave propagation, statistics of |srmt.2i| 
can be generated from RMT, and the only parameter of 
the distribution P{\srmt,2i\',j) is the dephasing rate 7 
defined in [IJ]. Hemmady et al. [ig] found the relation- 
ship between 7 and the loss parameter a of the corre- 
sponding closed system as 7 = Ana. The loss parameter 
a is the ratio of the closed-cavity mode resonance 3-dB 
bandwidth to the mean spacing between cavity modes, 
a = f/{2QAf). Here / is the frequency, Q is the typ- 
ical quality factor, and A/ is the mean spacing of the 
adjacent eigenfrequencies. For an open fading system, 
we consider an equivalent closed system in which uni- 
form absorption accounts for wave energy lost from the 
system, and we assume that we can define an equivalent 
loss parameter a for the open system [iGj. 

We can analytically derive the distribution of the fad- 
ing amplitude P{\srmt,2i\', ct) for special cases of a. For a 
lossless system (a = 0), the distribution of fading ampli- 
tude is a uniform distribution between and 1. For high 
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FIG. 1: The probability density functions P(|srmt,2i|) gen- 
erated from two theoretical models. Solid curves show the 
numerical results from the RMT model with different loss pa- 
rameters (a — 0, 0.1, 1, 10). For higher loss cases (a = 1 and 
a — 10), the corresponding Rayleigh distributions are shown 
as dashed curves. 



loss systems (a ^ 1), we can prove that the distribu- 
tion of fading amplitude goes to a Rayleigh distribution 

P{x — \srmt,2i\', c) = ^ exp ( -=^ j with the relationship 
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This result reveals the physical meaning of the a parame- 
ter of the Rayleigh fading model, which assumes that the 
real and imaginary parts of the complex quantity 5*21 are 
independent and identically distributed Gaussian vari- 
ables with mean and variance ct^. 

In Fig. [T] we illustrate P{\srmt,2i\',ce) for different 
loss parameter values from the RMT fading model. For 
higher loss cases, we also plot the corresponding Rayleigh 
distributions from Eq. ([Ij to show the convergence of the 
two models in the high loss limit. Note that distributions 
from the RMT model in the low loss region {a < 0.1) de- 
viate from a Rayleigh distribution. 

To apply wave chaos theory to practical systems, in 
addition to the universal Srmt.2i, we also need to account 
for the non-chaotic features of the wave system. We em- 
ploy the random coupling model (RCM) [la - flSj . which 
combines universal fluctuating properties of a scatter- 
ing system with the non-universal features arising from 
the port geometry and short-orbits between the ports, in 
the impedance matrix description. The impedance (Z) 
matrix specifies the linear relationship between the port 
voltages and the port currents, and is related to the scat- 
tering matrix hy Z = Zq(\ + S)/(l — S), where Zq is 
a diagonal matrix with elements equal to the character- 
istic impedances of the transmission lines connected to 
the ports. Another method to deal with the coupling be- 
tween the scattering channels and the system is the Pois- 
son kernel that presents the non-chaotic features as an 
average S in the scattering matrix description [13, [20] ■ 
The advantage of the RCM is that it can separate the 
chaotic (fluctuating) part and the non-chaotic (average) 
part in a simple additive format; 



2^{the) jg ^^ theoretical prediction of the raw measured 
impedance matrix Z, and Ravg and ^avg are the real and 
imaginary parts of the system-specific ensemble-averaged 
impedance matrix Zavg- The matrix Zavg can be ap- 
proximated by taking the average of the impedances of 
all realizations in a finite ensemble, Zavg = Z. The 
chaotic part is the universal impedance matrix z^mt = 

(1 -I- Srmt)/(1 - Srmt)- 

In the extended RCM [OjllMli Hart et al. express the 
system-specific features in the ensemble wave-scattering 
system as 
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where Rrad and X^ad are the real and imaginary parts of 
the radiation impedance matrix Zrad, which is a diagonal 
matrix that quantifies the radiation and near-field char- 
acteristics of the ports. The other system-specific feature 
is short (major) trajectory orbits. We define an orbit as a 
ray trajectory that originates from one port, bounces on 
the boundary of the system or on scattering objects, and 
then reaches a port. Note that the line-of-sight signal be- 
tween the two ports is the first (shortest) orbit, and that 
a short orbit is distinct from a periodic orbit in a closed 
system [2l|. We can compute the short-orbits contribu- 
tion matrix Zso of the Af shortest orbits from the known 
geometry in each realization of the ensemble [17, 18]. The 

matrix elements z[l[l^^ = X)* =i P^'^S exp(«'S'5^) is 
the sum of the short-orbit terms of the M shortest orbits 
between port a and port b. In each term we consider 
scattering on dispersing surfaces (in the geometry factor 



(m). 



Cai), the survival probability {P^'h) of the orbits due 
to the presence of mobile perturbers, and the propagation 

phase advance and loss (in the action S)^ ) . The contri- 
bution of trajectory orbits decreases exponentially with 
the orbit length 17, Hi- In a lossy system ZavJ -^ Zavg 
as M increases 13, 111], and only a limited number of 
short orbits are required to represent system-specific fea- 
tures that survive the ensemble average. 

According to RMT, the universal complex parameter 
Srmt,2i has zero mean, but the system-specific features 
of the sum of short orbits Zso brings about a non-zero 
bias in the impedance matrix Zavg (or Zavg)- Therefore, 
the measured S21 can have a non-zero mean, and this is 
similar in character to the Rice model. The Rice fad- 
ing model use the distribution P(x = ] s„nt. 21 ]; cr, i^) = 
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parameter z/ (i/ — > recovers the Rayleigh distribution), 
and /o(-) is the modified Bessel function of the first kind 
of order zero. The Rice model is an extension of the 
Rayleigh model in which the real and imaginary parts 
of 521 are still independent and identical Gaussian vari- 
ables with variance ct^, but the means are generalized to 
a biased mean of magnitude v. The Rice fading model 
is used in environments where one signal path, typically 
the line-of-sight signal, is much stronger than the others 
[i S Q 1 ^^d the 1/ parameter is related to the strength of 



the strong signal. More generally, we find that the RMT 
fading model in the high loss limit yields an explicit ex- 
pression for V in terms of the short orbit impedance 
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This result generalizes the meaning of v to include the in- 
fluence of all major (short) paths. Note that when there 
is a single strong signal that dominates the sum of all 
paths, the v parameter reverts to the original interpreta- 
tion of Rice fading. 

We have carried out experimental tests of the RMT 
fading model by measuring the complex 2x2 scattering 
matrix S in two quasi-two-dimensional ray-chaotic mi- 
crowave cavities. Both of these cavities have two coupling 
ports, which we treat as a transmitter and a receiver. Mi- 
crowaves are injected through each port antenna attached 
to a coaxial transmission line of characteristic impedance 
Zq = 50f2, and each antenna is inserted into the cavity 
through a small hole in the lid, similar to previous se- 
tups [la, E, [24I . The waves introduced are quasi-two- 
dimensional for frequencies below the cutoff frequency for 
higher order modes ('^ 19 GHz) due to the thin height of 
the cavities (8 mm in the z-direction). 

Classical ray-chaos arises from the shape of the cav- 
ity walls. One cavity is a symmetry-reduced "bow-tie 
billiard" made up of two straight walls and two circu- 
lar dispersing walls [la, [3 shown in Fig. [2jc), and the 
other cavity is a "cut-circle billiard" 22] shown in Fig. 
[2I^d). The scales of the billiards compared to the wave- 
lengths of the microwave signals (1.7 — 5.0 cm) put these 
systems into the semiclassical limit. To create an ensem- 
ble for statistical analysis, we add two cylindrical metal 
perturbers to the interior of the 1/4 bow-tie cavity and 
systematically move the perturbers to create 100 differ- 
ent realizations. For the cut-circle cavity, the perturber 
is a Teflon wedge that can be rotated inside the cavity. 
We rotate the wedge by 5 degrees each time and create 
72 different realizations. The perturbers can be consid- 
ered as scattering objects in the propagation medium, 
so changing the positions creates the equivalent of time- 
dependent scattering variations that give rise to fading. 
The 1/4 bow-tie cavity is made of copper, and measure- 
ments of the transmission spectrum at room temperature 
suggest the loss parameter goes from a = 0.4 to 1.0, vary- 
ing with the frequency range JJ,]. The superconducting 
cut-circle cavity is made of copper with Pb-plated walls 
and cooled by a pulsed tube refrigerator to a temper- 
ature of 5.5 K, below the transition temperature of Pb 
[23, [231 . Measurements of the transmission spectrum sug- 
gest a < 10~^. 

The experimental data show good agreement with our 
RMT fading model. We first use Z^*'*^) in Eq. © to 
represent the measured impedance matrix Z and solve 
for Zrmt (and therefore Srmt)- In this process we re- 
move the system-specific features including all short or- 
bits, so the situation is equivalent to the Rayleigh fad- 
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FIG. 2: Probability density functions P{\srmt,2i\) from the 
experimental data (red circles) in (a) the 1/4 bow-tie cav- 
ity and in (b) the cut-circle cavity, comparing with the RMT 
model (black solid) and the best-matched Rayleigh distribu- 
tion (blue dashed), (c) The 1/4 bow-tie cavity with the two 
ports as red dots and the two perturbers as blue circles, (d) 
The cut-circle cavity with the two ports as red dots and the 
perturber as a blue wedge, (e) Magnitude of |sl,o 21I averaged 
over a 2-GHz frequency band versus the central frequency 
in the 1/4 bow-tie cavity, and the u parameter of the best- 
matched Rice distribution. 



ing environment where no direct paths exist. By choos- 
ing data over all realizations in a frequency range, we 
can construct the distribution of |srmt.2i| and compare 
with the prediction of RMT. In Fig. [1 (a) and (b) we 
plot the distributions of the fading amplitude from the 
RMT model (black solid), the experimental data (red 
circles), and a best-matched Rayleigh distribution (blue 
dashed). In Fig. ^a) the room-temperature case, the 
best-matched RMT model gives a value of the loss pa- 
rameter a — 0.5 for the experimental data, which cor- 
responds to cr ~ (87ra)~"-^ = 0.282. The best-matched 
Rayleigh distribution yields a = 0.226. The difference 
in a values is due to the fact that the loss parameter 
is not very large in this case. Nevertheless, both mod- 
els agree with the experimental data well in this loss 
regime. In Fig. ^h) the superconducting cavity case, 
the agreement between the experimental data and the 
RMT model is much better than the Rayleigh distribu- 
tion. In fact, in the very-low- loss region (a << 1), the 
long exponential tail of a Rayleigh distribution can never 
match the RMT theoretical distribution that is limited 

to < \s„nt,2l\ < 1- 

In the room-temperature case, since the loss parameter 
is high enough, we can compare the relationship between 
the RMT model and the Rice fading model [Eq. (g])]. 
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In Fig. mje), we compute 

with length up to 200 cm (M = 1088) in the 1/4 bow-tie 



cavity, apply Eq. (j4]), perform a sliding average over a 

2-GHz frequency band, and plot the magnitude |Sso2il 
as the red curve. For the v parameter of the Rice model, 
we first remove the coupling features from the measured 

impedance (Z) matrix as z = R~^^ {Z - iXrad)R~ad 
and convert the impedance matrix z to s. Then we com- 
pare the distribution of |s2i | over a 2-GHz frequency band 
and 100 realizations with the best-matched Rice distribu- 
tion. Since the a parameter has been determined by the 
best-matched Rayleigh distribution as described above 
for the fully universal data [Fig. [D^a)], we can use ly as 
the only fitting parameter. In Fig. ^e) we plot the i/ 
parameter for the best-matched Rice distributions (blue 
squares) along with the system-specific average magni- 
tudes of S21 versus the central frequency of a 2-GHz fre- 
quency band. The value of the Rice v parameter and 
the system-specific feature described by our model agree 
well. 

One more advantage of applying the RCM is that we 
can extend the relations Eq. ^ and (|4]) from the normal- 
ized data to the raw measured data in the high loss cases. 
In high loss cases, the magnitude of the elements of Srmt 
are much less than one, so we take the approximation to 
the lowest order [2^. For the generalized D parameter, 

we only need to replace the Zso terms by Zavg or Zavg 
in Eq. Q. The generalized a parameter is a function of 



a and all elements of the matrix Zavg ■ If the transmission 
between the ports is much less than the coupling reflec- 
tion at the ports (i.e. \Zavg,2l\ < \Zavg,ll\, \Zavg,22\), 

the modified Rayleigh a parameter can be simplified to 
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In conclusion, we have provided a first-principles 
derivation of a RMT fading model that reduces to the 
traditional Rayleigh and Rice fading models in high-loss 
scattering environments, and hence we can explain the 
physical meanings of the a parameter of the Rayleigh dis- 
tribution and the v parameter of the Rice distribution. 
Moreover, in low-loss environments, the RMT model can 
better predict the distribution of the fading amplitude 
5*21!. Because wave propagation in a complicated envi- 
ronment is a common issue in many fields [l|-|5|, |7[ , the 
fading model can be applied to wireless communication, 
global positioning systems, and mesoscopic physics. 
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